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domain of human coronavirus spike 
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Continuously emerging highly pathogenic human coronaviruses (HCoVs) remain a major threat to human health, 
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as illustrated in past SARS-CoV and MERS-CoV outbreaks. The development of a drug with broad-spectrum HCoV 
inhibitory activity would address this urgent unmet medical need. Although previous studies have suggested that 
the HR1 of HCoV spike (S) protein is an important target site for inhibition against specific HCoVs, whether this 
conserved region could serve as a target for the development of broad-spectrum pan-CoV inhibitor remains con- 
troversial. Here, we found that peptide OC43-HR2P, derived from the HR2 domain of HCoV-OC43, exhibited broad 
fusion inhibitory activity against multiple HCoVs. EK1, the optimized form of OC43-HR2P, showed substantially im- 
proved pan-CoV fusion inhibitory activity and pharmaceutical properties. Crystal structures indicated that EK1 can 
form a stable six-helix bundle structure with both short a@-HCoV and long B-HCoV HR's, further supporting the 


role of HR1 region as a viable pan-CoV target site. 


INTRODUCTION 

Coronaviruses (CoVs) are enveloped viruses with a positive-sense, 
single-stranded RNA and are associated with various natural hosts. 
CoVs are divided into alpha, beta, gamma, and delta groups, and 
the beta group is further composed of A, B, C, and D subgroups 
(Fig. LA) (1). Among them, six CoVs can infect humans (HCoVs), 
including HCoV-229E (229E) and HCoV-NL63 (NL63) in the alpha 
group, HCoV-OC43 (OC43) and HCoV-HKU1 (HKU1) in beta 
subgroup A, severe acute respiratory syndrome CoV (SARS-CoV) 
in beta subgroup B, and Middle East respiratory syndrome CoV 
(MERS-CoV) in beta subgroup C (2). 

In this century, SARS-CoV and MERS-CoV have emerged in the 
human population and caused severe pulmonary disease with 
alarmingly high case-fatality rates. In 2002, SARS-CoV infections 
first appeared in China and then quickly spread as a global epidemic 
in more than 30 countries with 8273 infections and 775 deaths (nearly 
10% mortality) (2). In 2012, MERS-CoV emerged in Saudi Arabia 
and spread throughout the Middle East. In 2015, the second pandemic 
of MERS-CoV occurred in South Korea, causing super-spreading 
events with third- and fourth-generation cases of infection. The 
World Health Organization has reported 2229 laboratory-confirmed 
cases of MERS-CoV infection, including 791 deaths (about 35% case 
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fatality) in 27 countries as of August 2018 (www.who.int/emergencies/ 
mers-cov/en/). Meanwhile, the remaining common HCoVs, such as 
229E, OC43, and NL63, usually infect the human upper respiratory 
tract and cause the common cold, but they also are responsible for 
severe and even fatal diseases in children, elderly, and immuno- 
compromised patients (3-5). These scenarios suggest that those 
common HCoVs might also pose a lethal threat to humans. Note 
that HCoVs are rapidly evolving. OC43 isolates with novel genomes 
are being continuously identified (6-8). 

The existence of SARS-like CoV (SL-CoV) and MERS-like CoV 
(ML-CoV) also pose great threats to public health worldwide. Re- 
cent studies identified some types of SL-CoV, such as SL-WIV1-CoV 
and SL-SHC014-CoV, whose spike (S) proteins highly resemble those 
of SARS-CoV. These SL-CoVs could use the same SARS-CoV recep- 
tor, 1.e., angiotensin-converting enzyme-2 (ACE2), to directly enter 
permissive human cells without need for adaptation (9-11). In addi- 
tion, the ML-CoV, Tylonycteris bat CoV HKU4, was shown to recog- 
nize the MERS-CoV receptor CD26 and infect human cells either after 
introduction of two single mutations (S746R and N762A) into its S 
protein or with the help of exogenous protease (12, 13). From a histor- 
ical perspective, zoonotic CoVs harbor a strong ability to cross species 
barriers to infect humans rapidly and unpredictably, as illustrated by 
newly emerging HCoVs (2, 14). Thus, developing a specific drug that 
only targets a single HCoV would be ineffective against newly emerg- 
ing HCoVs (9-11). Since no broad-spectrum anti- HCoV drug is cur- 
rently available for clinical use, it is incumbent to search for a common 
or conserved target site based on existing HCoVs. 

The S glycoprotein is a type I transmembrane glycoprotein that 
plays an important role in mediating viral infection and is common 
to all HCoVs. The S proteins consist of two subunits, S1 and S2 
(Fig. 1B). The S1 subunit binds the cellular receptor through its 
receptor-binding domain (RBD), followed by conformational changes 
in the $2 subunit, which allows the fusion peptide to insert into the 
host target cell membrane. The heptad repeat 1 (HR1) region in 
the $2 subunit forms a homotrimeric assembly, which exposes 
three highly conserved hydrophobic grooves on the surface that 
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Fig. 1. HR1P and HR2P fusion inhibitory activity. (A) Maximum likelihood trees created with the S sequence of representative CoVs from all four genogroups. HCoVs 
are denoted with filled red circles. (B) The antiviral mechanism of HR2P peptides. (C) Schematic representation of HCoV S protein. SP, signal peptide; FP, fusion peptide; 
HR, heptad repeat domain (HR1 and HR2); TM, transmembrane domain; CP, cytoplasmic domain. Corresponding sequences of the designed peptides (HR1Ps and HR2Ps). 
(D) HR2P peptides potently inhibit cell-cell fusion mediated by the S proteins of MERS-CoV, SARS-CoV, 229E, NL63, and OC43. Data are means + SD of triplicate samples from 


a representative experiment. 


bind heptad repeat 2 (HR2). This six-helix bundle (6-HB) core structure 
is formed during the fusion process and helps bring the viral and 
cellular membranes into close proximity for viral fusion and entry 
(Fig. 1B) (15). Thus, the S protein is an important target protein for 
the development of specific drugs. In particular, the S1 RBD is a very 
good target site, and both RBD-specific antibodies and RBD-based 
vaccines have previously exhibited effective antiviral activity or pro- 
tective effect in blocking binding of virus to host receptors (16). How- 
ever, from an evolutionary perspective, the RBD of CoV is part of a 
highly mutable region and, thus, is not an ideal target site for broad- 
spectrum antiviral inhibitor development (14). The SARS-CoV RBD- 
specific antibody fm6 failed to block infection mediated by the S protein 
of SL-CoV-SHC014 (9-11, 17). In contrast, the HR region in the 
S2 subunit is conserved among various HCoVs and plays a pivotal 
role in HCoV infections by forming the 6-HB that mediates viral 
fusion (fig. S1A). Furthermore, the mode of interaction between 
HR1 and HR2 is conserved among CoVs such that residues located 
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at the “e” and “g” positions in the HRI helices interact with residues 
at the “a” and “d” positions in the HR2 helices (fig. S1B) (18). Previ- 
ous studies have reported that peptides derived from the HR2 (or 
C-terminal HR) region of class I viral fusion proteins from some en- 
veloped viruses, including HIV-1 (19-21), respiratory syncytial virus 
(RSV) (22), Ebola virus (23), paramyxoviruses SV5 (24), Nipah 
virus (25), and murine hepatitis virus (MHV) (26), could com- 
petitively bind the viral HR1 (or N-terminal HR) and effectively in- 
hibit viral infection. Therefore, it is reasonable to speculate that HR1 
could also be a good target for the development of fusion inhibitors 
against highly pathogenic HCoVs. 

We and others have reported that peptides derived from the 
HR2 regions of SARS-CoV and MERS-CoV S proteins can competi- 
tively inhibit viral 6-HB formation, thereby preventing viral fusion and 
entry into host cells (18, 27). For example, the peptide CP-1, derived 
from the SARS-CoV spike protein HR2 region, was able to in- 
hibit SARS-CoV entry in a manner similar to that of MERS-HR2P 
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against MERS-CoV infection (18, 27). However, we note that those 
peptides lack broad-spectrum antiviral activity against heterolo- 
gous HCoVs. For example, CP-1 and MERS-HR2P peptides failed 
to cross-inhibit MERS-CoV and SARS-CoV infection, respectively. 
Furthermore, the CoV fusion core HR regions can be divided 
into two groups (28, 29): short HRs, such as MERS-CoV, SARS- 
CoV, and OC43 HRs; and long HRs, such as 229E and NL63 HRs. 
The difference between short and long HRs arises from an in- 
sertion of 14 amino acids in the long HRs, further increasing the 
difficulty in designing broad-spectrum peptide fusion inhibitors. 
To address this challenge, we herein report the successful screening 
of a peptide OC43-HR2P with broad-spectrum fusion inhibitory 
activity. Furthermore, a modified OC43-HR2P peptide, EK1, shows 
promising potency and breadth in inhibiting infection by multi- 
ple HCoVs. In vivo studies demonstrate that administration of 
EK1 via the nasal route exhibits highly protective effects and 
safety profiles, highlighting its clinical potential. Moreover, 
structural studies of EK1 in complex with HR1s from different 
HCoVs explain the conserved basis for the HR1-EK1 interaction, 
further indicating that HR1 region could serve as a promising tar- 
get site for the development of broad-spectrum pan-CoV fusion 
inhibitors. 


RESULTS 

HR-derived peptides inhibit cell-cell fusion mediated by S 
proteins of multiple HCoVs 

On the basis of the 6-HB fusion core structure, we previously re- 
ported on peptides MERS-HRI1P and MERS-HR2P that were derived 
from MERS-CoV HRI and HRz2, respectively (fig. S1C), and these 
two peptides displayed good interaction with each other (27). Here, 
we located the conserved HR region of multiple HCoVs by sequence 
comparisons and then synthesized HR1- and HR2-derived peptides, 
termed HR1P and HR2P (Fig. 1C). Notably, we observed 14-amino 
acid insertions in both HR1P and HR2P in the HRs of two o-HCoVs, 
i.e., 229E and NL63. 

To systematically assess the inhibitory activities of these peptides 
against different HCoVs, we developed multiple cell-cell fusion 
assays that were mediated by the S protein of various HCoVs (fig. 
S1D). Consistent with previous results, MERS-HR2P had inhibitory 
activity against cell-cell fusion mediated by MERS-CoV S protein 
with a concentration for 50% inhibition (ICs9) of 1.01 uM, whereas 
it showed little inhibitory activity in other HCoV S-mediated cell- 
cell fusion assays even with concentrations up to 5 uM (Fig. 1D and 
table S1). Similarly, SARS-HR2P specifically inhibited SARS-CoV 
S-mediated cell-cell fusion with an ICs9 of 0.52 uM. On the other 
hand, HR2P peptides derived from the two a-HCoVs, i.e., 229E and 
NL63, showed potent and broad inhibitory activity against a-HCoV 
S-mediated cell-cell fusion with ICso values ranging from 0.13 to 
0.51 uM or from 0.21 to 0.56 uM, respectively, but no effective 
inhibitory activity against B- HCoVs (including MERS-CoV, 
SARS-CoV, and OC43) S-mediated fusion (Fig. 1D and table S1). 
OC43-HR2P exhibited broad and potent fusion inhibitory activity 
with ICs values of 0.39, 0.54, and 0.66 uM against MERS-CoV, 
SARS-CoV, and OC43, respectively (Fig. 1D and table S1). Un- 
expectedly, OC43-HR2P, a B-HCoV HR2-derived peptide, and thus 
14 residues shorter than 229E-HR2P and NL63-HR2P, exhibited 
effective activity against o-HCoVs with an ICs 9 of 0.84 uM on 
229E-S—mediated cell-cell fusion and an ICs9 of 0.94 uM on NL63-S- 
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mediated cell-cell fusion. Among all HR1Ps, only 229E-HRIP ex- 
hibited moderate inhibitory effects on 229E-S— and NL63-S—mediated 
cell-cell fusion (Fig. 1D). Thus, compared to HRIPs, all HCoV 
HR2-derived peptides exhibited excellent self-specific fusion in- 
hibitory activity, whereas OC43-HR2P showed broad-spectrum 
and potent fusion inhibitory activity against both o-HCoV and 
B-HCoV S-mediated cell-cell fusion (Fig. 1D and table S1). We also 
measured the a-helicity of OC43-HR2P and HCoV-HRIPs by circu- 
lar dichroism and found that the peptides alone exhibited limited 
a-helicity ranging from 11.3 to 42.3% (fig. SIE). In contrast, the 
mixtures of OC43-HR2P and each of these HR1 peptides, respec- 
tively, exhibited high a-helicity (70.9 to 86.7%) with melting tran- 
sition temperature (T,,) values that ranged from 48.5° to 91.5°C 
(fig. S1E), further implying that OC43-HR2P can bind HRIs of 
different HCoVs and form stable complexes, thereby blocking 
S protein-mediated fusion. 


EK1, a modified OC43-HR2P peptide, has improved fusion 
inhibitory activity and increased solubility 

It was previously reported that introduction of negatively and posi- 
tively charged amino acids Glu (E) and Lys (K), at the itoi+3 ori 
to i+ 4 positions in a helix, into an HIV-1 fusion inhibitory peptide 
can form intramolecular E-K or K-E salt bridges that result in sub- 
stantial enhancement of the peptide’s stability, solubility, and anti- 
viral activity (30, 31). Using a similar design, we further optimized 
the sequence of OC43-HR2P by introducing Glu or Lys at appropriate 
sites in the peptide to increase the solubility and thereby the antiviral 
activity of the peptide. Moreover, on the basis of the structure of 
MERS-CoV S 6-HB, we also introduced mutations at some sites that 
are not expected to be involved in HR1 binding, such as 4Q, 14Y, 
32D, and 36L, to further enhance the fusion inhibitory activity of 
the peptide (table S2). Among the series of optimized peptides, pep- 
tides EKO-1, EKO-2, EKO-3, and EK1 showed gradually increased 
solubility and excellent inhibitory activity in cell-cell fusion assays. 
The final peptide, EK1, exhibited the most potent pan-CoV antiviral 
fusion activity with IC59 values in the range of 0.19 to 0.62 uM (table 
$2). In some HCoV cell-cell fusion assays, EK1 exhibited even more 
effective inhibitory activity than the autologous peptide, such as 
MERS-HR2P and SARS-HR2P. Akin to its ancestor OC43-HR2P, 
EK1, which is derived from the short 6-HB fusion core of OC43, was 
also able to potently inhibit the cell-cell fusion mediated by 229E and 
NL63 S proteins, both of which harbor the long fusion core. These 
results further underscore the broad-spectrum anti-HCoV potential 
of EK1 (Fig. 2, A to E, and table $2). Furthermore, EK1 also had 
superior pharmaceutical properties and solubility in phosphate- 
buffered saline (PBS) and water, which increased by 478-fold and 
3.5-fold, respectively, as compared to the original OC43-HR2P. 
Together, these results establish an important framework and plat- 
form for subsequent development of EK1 as a potential therapeutic 
(table S2). 


EK1 peptide potently inhibits multiple CoV cell-cell fusion 
and blocks various pseudotyped and live CoV infection 
SL-CoVs, including WIV 1, Rs3367, and SHC014 CoVs, all manifest 
potential for human infection (9-11). To further assess the breadth 
of fusion inhibitory activity, as demonstrated by EK1, we estab- 
lished cell-cell fusion assays mediated by the S protein of these 
SL-CoVs (fig. S2). Notably, many studies have suggested diversity in 
the receptor-binding motif (RBM) of SL-CoV’s RBD (9-11), which 


3 of 15 


6LOZ ‘OL |udy uo /Bio bewseousios ssoueApe//:dy1y Wo) pepeojuMOG 


SCIENCE ADVANCES | RESEARCH ARTICLE 


CoVs S-mediated cell-cell fusion 


A MERS-Cov DB SARS-CoV Cc HCov-229E D Hcov-ntes SClCi HCov-oc43.—O#F Rs3367 G wiv H SHC014 
== EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled 
== MERS-HR2P == SARS-HR2P == 229E-HR2P = NL63-HR2P == 0C43-HR2P = SARS-HR2P == SARS-HR2P == SARS-HR2P 
1004 EK 1004 “™ EK 1ooq “™ EK1 1004 “™ EK1 1004 ™ EK1 _ 100 = 1004 EK tog EK 
S 0 S80 = 80 = 80 = 80 S 80 80 80 
5 60 5 60 & 6 5 60 & 60 6 60 5 6 5 60 
5 40 6 40 6 40 5 40 5 40 5 40 5 40 6 40 
= 2 = 2 = 2 = 2 = 2 = 2 = 2» = 2 
0 0 0 0 0 0 0 0 
0.1 4 10 0.1 1 10 0.1 1 10 0.1 1 10 0.1 4 10 0.1 1 10 0.1 1 10 0.1 1 10 
Concentration (uM) Concentration (uM) Concentration (uM) Concentration (uM) Concentration (uM) Concentration (uM) Concentration (uM) Concentration (uM) 
l MERS-cov J SARS-CoV K HCov-229—E 3=6ssL Hcov-nL6es3. M Hcov-ocaa. ~=ON Rs3367 O Wwiv1 P VSV 
== EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled == EK1-scrambled ~~ EK1 
=~ MERS-HR2P == SARS-HR2P =® 229E-HR2P = NL63-HR2P = 0C43-HR2P =~ SARS-HR2P == SARS-HR2P «= EK1-scrambled 
1004 EK 100g ™ EK! 100g ~~ EK oo. 1004 ™ EK 1004 EKA we 400 
= a0 = 80 S 80 ss 80 Ss 80 = 80 S 30 S 0 
& 60 & 60 & 60 § 60 § 60 6 60 & 60 & 60 
5 40 5 40 5 40 5 40 5 40 5 40 6 40 S40 
= 20 = 20 = 2 = 2 = 2 = 2 = 2 = 2 
0 0 0 0 0 0 0 0 F : i 
0.1 1 10 0.1 1 10 0.1 1 10 0.1 1 10 0.1 1 10 0.1 1 40 0.1 1 10 0.4 4 40 
Concentration (uM) Concentration (uM) Concentration (uM) Concentration (1M) Concentration (uM) Concentration (uM) Concentration (uM) Concentration (uM) 
Q MERS-CoV R HCoV-0C43 S HCoV-229E T HCoV-NL63 
E 5 ea cka E 6 @m EK1 = qm EK1 £ @m EK1 
3 ro @™ 0C43-HR2P 3 229E-HR2P B @™ NL63-HR2P 
Q @ MERS-HR2PM2 a: xg a oO g 4 
O oO oO Ss) 
= E 
. ° 4 e* o 3 
D Do 
° re) fe) 
= = 3 = — 2 
= — 
2 2 5 £ 2 
+ boas ear ~ 1 
= 5 _ i 
£ = 1 a £ 
> 0.00 0.16 0.63 2.60 10.00 20.00 a 0.00 0.16 0.63 2.50 10.00 20.00 = 0 016 063 25 10 #20 > 0 0.02 0.08 0.31 1.25 5 20 


Concentration (uM) Concentration (uM) 


Concentration (uM) 


Concentration (uM) 


Fig. 2. EK1 is effective against viral infection mediated by S protein of multiple HCoVs. (A to H) Inhibitory activity of EK1 in cell-cell fusion mediated by the S proteins 
of MERS-CoV (A), SARS-CoV (B), 229E (C), NL63 (D), OC43 (E), Rs3367 (F), WIV1 (G), and SHCO14 (H). (I to P) Inhibitory activity of EK1 in pseudovirus infection assays against 
MERS-CoV (I), SARS-CoV (J), 229E (K), NL63 (L), OC43 (M), Rs3367 (N), WIV1 (O), and VSV (P). (Q to T) Inhibitory activity of EK1 on live HCoV replication for MERS-CoV (Q), 
OC43 (R), 229E (S), and NL63 (T). Data are means + SD of triplicate samples from a representative experiment. 


plays an important role in binding the cellular receptor and serves 
as an effective target site for the development of CoV-specific anti- 
bodies or vaccines (32). Specifically, the RBM of SHC014 CoV has 
only 80% similarity and 64% identity to that of SARS-CoV, and an 
antibody targeting the SARS-CoV RBD could not prevent infection 
mediated by the SHC014 S protein (9-11, 17). In contrast, the HR1 
and HR2 sequences of SL-CoV are 100% identical to those of SARS- 
CoV (fig. S2B). Correspondingly, EK1 exhibited greater inhibitory 
activity than the autologous peptide SARS-HR2P against cell-cell 
fusion mediated by the S protein of the three SL-CoVs tested 
(Fig. 2, F to H), while an EK1-scrambled peptide showed no inhibi- 
tory activity on cell-cell fusion mediated by any of these CoV S 
proteins (Fig. 2, A to H). 

The pseudovirus assay is a good model to mimic the process of 
virus entry into the target cell and has been widely used in previous 
studies to assess the inhibitory activity of antiviral agents against 
related CoV infection (33). We also used here pseudovirus assays to 
assess the inhibitory activity of EK1 against different pseudotyped 
HCoVs. Ina MERS-CoV pseudovirus infection assay, both EK1 and 
MERS-HR2P showed inhibitory activity with ICso values of 0.26 
and 1.06 uM, respectively (Fig. 21). Similarly, EK1 and SARS-HR2P 
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also effectively blocked SARS-CoV pseudovirus infection with ICs 
values of 2.23 and 2.81 uM, respectively (Fig. 2J). Meanwhile, EK1 
demonstrated antiviral activity against pseudotyped 229E, NL63, 
and OC43 infection with ICs values of 3.35, 6.02, and 1.81 uM, re- 
spectively (Fig. 2, K to M). Consistent with the results from cell-cell 
fusion assays, both EK1 and SARS-HR2P inhibited pseudotyped 
Rs3367 virus infection with ICs9 values of 2.25 and 3.05 uM, respec- 
tively, and prevented WIV1 pseudovirus infection with ICs9 values 
of 2.10 and 2.73 uM, respectively (Fig. 2, N and O). In contrast, EK1 
exhibited no inhibitory activity against VSV-G-mediated viral 
infection (Fig. 2P). The EK1-scrambled peptide showed no 
inhibitory activity on infection by any of these pseudoviruses 
(Fig. 2, I to P). As previously reported (10), we failed to assemble 
a pseudotyped SHC014 and thus did not perform pseudovirus in- 
fection assays on this virus. Collectively, pseudovirus infection by 
all HCoVs and SL-CoVs could be efficiently blocked by EK1, fur- 
ther indicating that EK1 has broad-spectrum antiviral activity 
against infection by pan-HCoVs, including a-HCoVs, B-HCoVs, 
and SL-CoVs. 

Blam-Vpr assay is a sensitive model to characterize the inhibitory 
mechanism and activity of antiviral agents on virus—cell membrane 
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fusion. This assay has been widely used for assessing the membrane 
fusion activity of various viruses, including HIV-1, Ebola virus, and 
henipavirus (34-37). Here, we developed the Blam-Vpr assays for 
MERS-CoV and SARS-CoV to further clarify the fusion inhibitory 
mechanism of EK1. We found that EK1 effectively inhibited the 
entry of MERS-CoV Blam-Vpr virions and SARS-CoV Blam-Vpr 
virions into Huh-7 cells and ACE2-293 T cells, respectively, in a 
dose-dependent manner (fig. S2D), thus confirming that EK1 does 
act as an HCoV fusion inhibitor. 

Subsequently, we further evaluated the inhibitory activity of EK1 
against live HCoV infection, including MERS-CoV, OC43, 229E, 
and NL63. We found that EK1 could effectively inhibit the infection 
and replication of these HCoVs at cellular level in a dose-dependent 
manner (Fig. 2, Q to T). EK1 even surpassed the autologous OC43- 
HR2P (ICs = 0.93 uM) or MERS-HR2M2 [ICs9 = 0.23 UM; which is 
an optimized MERS-HR2P peptide (33)] against OC43 or MERS- 
CoV with ICso values of 0.62 and 0.11 uM, respectively. Mean- 
while, EK1 also exhibited equally effective inhibitory activity 
against 229E and NL63 infections with ICso values of 0.69 and 
0.48 uM, respectively, compared to 229E-HR2P (ICs = 0.56 uM) 
and NL63-HR2P (ICs = 0.14 uM) (Fig. 2, S and T). Unfortunately, 
because of the global prohibition of handling live SARS-CoV in 
most virology laboratories, we could not test the inhibitory activity 
of EK1 against SARS-CoV replication. 


EK1 exhibits broad anti-HCoV activity in vivo 
The respiratory tract is the primary target tissue of HCoV infection. 
To investigate the preclinical potential of the EK1 peptide, we 
administered it via an intranasal route and assessed its distribution 
within the respiratory tract. We observed the fluorescent signals of 
Cy5-EK1 mainly in the upper and lower respiratory tract of mice 
(Fig. 3, A and B). Compared with those from PBS-treated mice, the 
lungs obtained from Cy5-EK1-treated mice exhibited significantly 
higher fluorescence signals (P < 0.0001, n = 3) (Fig. 3, Cand D). To- 
gether, these results suggest that EK1 (intranasal) can widely be dis- 
tributed in the whole respiratory tract and be enriched in the lung. 
We also noted that Cy5-EK1 administered intranasally could be de- 
tected in significant amounts in several extrapulmonary organs, includ- 
ing the liver, kidney, and spleen, in some animals, suggesting that EK1 
could enter into the blood circulation and other organs (fig. $3, A and 
B). This observation suggests that intranasal administration of EK1 
could also be beneficial for multiorgan infection or systemic infection 
of HCoVs, both of which are common in MERS-CoV infection. 
Next, we assessed the protective effect of EK1 in vivo on OC43 
and MERS-CoV infection mouse models. In the OC43-infected mouse 
model, we treated newborn mice with EK1 at a dose of 5 mg/kg or 
with PBS 30 min before or after challenge with HCoV-OC43 at 107 
TCIDs9 (50% tissue culture infectious dose). Body weight of mice in 
the viral control group decreased, starting from 5 days postinfection 
(dpi), and mice succumbed to infection by 10 dpi with 100% mor- 
tality (Fig. 3, E to G). In contrast, the final survival rate of mice in 
the EK1 prophylactic and therapeutic groups was 100 and 66.7%, 
respectively (Fig. 3E), and their body weight either appeared normal 
(prophylactic) or rapidly recovered at 16 dpi (therapeutic) (Fig. 3F). 
Meanwhile, we tested the viral titers in brains of mice of all groups 
at 5 dpi. Infectious virus was readily detectable in the viral control 
group, whereas infectious virus titers were below the limit of detec- 
tion (2 log TCIDs/g) in the EK1 prophylactic mice or very low in 
the EK1 therapeutic mice (Fig. 3G). However, the infectious virus 
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titers in the brains of mice that died during EK1 therapeutic treat- 
ment were as high as those in the brains of viral control mice with- 
out EK1 treatment, while the viral titers in the brains of survival mice 
with EK1 therapeutic treatment and in those of normal control 
mice were both undetectable (fig. S2E). Consistently, the brains of 
mice that died during EK1 therapeutic treatment exhibited similar 
histopathological changes as those of the viral control mice, i.e., simi- 
lar amount of vacuolation, degeneration, and infiltration (fig. S2F). 
In contrast, the brains of the survival mice with EK1 therapeutic 
treatment and those of the normal control mice showed no apparent 
histopathological changes (fig. S2F). 

To evaluate the prophylactic and therapeutic potentials of EK1 
against MERS-CoV infection, we took advantage of the well-characterized 
transgenic (Tg) mice globally expressing human dipeptidyl peptidase 
IV (DPP4) viral receptor. We treated mice with 200 ug of EK1 or 
with PBS 30 min before or after challenge with MERS-CoV at 10° 
TCIDs9. Likewise, the weight loss of Tg mice treated with EK1 before 
or after MERS-CoV challenge was insignificant and rapidly recov- 
ered at 16 dpi, while untreated mice progressively lost significant 
weight before succumbing to infection within 9 dpi with 100% 
mortality (Fig. 3, H to J). In contrast, the survival rates in the EK1 
prophylactic and therapeutic groups were 100 and 75%, respectively 
(Fig. 3, H to J). When the yield of infectious viruses in lungs was 
used as the end point for assessing the efficacy of EK-1 on MERS- 
CoV infection, we were unable to recover any infectious virus from 
both EK1 prophylactic and therapeutic groups, in sharp contrast to 
the untreated controls (Fig. 3J). Together, these results indicated 
that the peptide EK1 has broad and potent prophylactic and thera- 
peutic efficacy against HCoV infections. 


The safety and low immunogenicity of EK1 

Safety is obviously very important for the development of EK1 in 
clinical applications, and therefore, we first tested its in vitro cytotox- 
icity on various target cells. EK1 is not cytotoxic at concentrations 
up to 1 mM, which is more than 200 times higher than its ICso for 
inhibiting any HCoV S-mediated cell-cell fusion and pseudovirus 
entry (fig. S3C). We then further investigated its cytotoxicity in vivo. We 
continuously administered mice with PBS, low-dose EK1 (20 mg/kg), 
or high-dose EK1 (100 mg/kg) by intranasal route every day for 1 week, 
and we recorded the body weight changes for the following 2 weeks 
(fig. S3D). The EK1-treated mice in both high- and low-dose groups 
lived normally, with no apparent difference in weight gain/loss ob- 
served as compared to PBS-treated mice. Meanwhile, we used enzyme- 
linked immunosorbent assay (ELISA) to measure EK1-specific 
antibodies in sera of mice 2 weeks after intranasal administra- 
tion. No EK1-specific antibodies were detected in either the EK1- 
treated group or in the PBS-treated group (fig. S3E). 

On the basis of the distribution of EK1 by intranasal administra- 
tion (fig. S3, A and B), we postulated that EK1 could penetrate the 
air-blood barrier to enter blood circulation and become enriched in 
some important organs, such as the liver and kidney. Therefore, we 
further examined levels of alanine aminotransferase (ALT; fig. S3F) 
and creatinine (fig. S3G) in the sera of mice. ALT and creatinine 
showed no significant difference (P > 0.05) at all time points be- 
tween EK1- and PBS-treated groups, suggesting that nasal applica- 
tion of EK1 at high or low doses did not affect mouse hepatic and 
renal function. 

We further compared the potential histopathological changes 
of EK1- and PBS-treated groups 4 weeks after administration. The 
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Fig. 3. In vivo prophylactic and therapeutic efficacy of EK1 in mice against OC43 and MERS-CoV infection. (A and B) Imaging of mice treated with Cy5-EK1 or PBS 
by the IVIS Lumina K Series III from PerkinElmer and the statistical analysis. (C and D) Imaging of lungs from those mice with the statistical analysis. (E to G) Anti-OC43 
efficacy of EK1 in vivo. Survival curves of 3-day-old suckling mice challenged with OC43. (E) Newborn mice were treated with EK1 (5 mg/kg in PBS) or PBS 30 min before 
or after challenge with OC43 (10* TCIDs9). (F) Body weight change of newborn mice treated with EK1 or PBS 30 min before and after OC43 challenge. (G) Viral titer in 
mouse brain of each group. (H to J) Anti-MERS-CoV efficacy of EK1 in vivo. Survival curves of mice challenged with MERS-CoV. (H) Mice expressing human DPP4 were 
treated with 200 wg EK1 in PBS or PBS 30 min before or after challenge with MERS-CoV (107 TCIDs0). (I) Body weight change of each group of mice. (J) Viral titer in mouse 
lungs of each group. SD of triplicate samples from a representative experiment; *P < 0.05, ***P < 0.001, and ****P < 0.0001. 


hematoxylin and eosin staining of lung, liver, kidney, and spleen 
sections from mice treated with EK1 at different doses exhibited no 
pathological abnormality, when compared to control mice treated 
with PBS (fig. S3H). None of these organs showed evidence of cell 
degeneration, necrosis, or infiltration of inflammatory factors. Hence, 
EK1 appears to be generally safe via nasal application. 


Structural basis of EK1 pan-CoV inhibitory activity 

To investigate the structural basis for the pan-CoV inhibitory effect 
of EK1 peptide, we crystallized EK1 in complex with HR1 peptides 
from three representative HCoVs, including the most pathogenic 
SARS-CoV and MERS-CoV in B-HCoVs and 229E in a-HCoVs, the 
HRI of which is 14 amino acids longer than those of B-HCoVs. 


Xia et al., Sci. Adv. 2019; 5:eaav4580 10 April 2019 


Crystal structures of HR1(MERS)-L6-EK1, HR1(SARS)-L6-EK1, 
and HR1(229E)-L6-EK1 were solved by molecular replacement and 
refined to final resolutions of 3.3, 3.3, and 2.2 A, respectively (table S3). 
In all three structures, EK1 snugly fits into the hydrophobic groove 
formed between two adjacent HR1 helices in an oblique and antiparal- 
lel manner (Fig. 4, A to C), producing a 6-HB structure resembling 
the 3HR1-3HR2 postfusion state of corresponding HCoVs (fig. S4). 
Notably, the binding sites of EK1 coincide with those of native HR2s 
(fig. $4). These results indicate that the presence of EK1 would preclude 
binding of HR2 onto their corresponding 3HRI core, thereby blocking 
formation of the 3HR1-3HR2 6-HB, which is an indispensable step 
during host-viral membrane fusion. Hence, administration of the EK1 
peptide would block cellular entry of those HCoVs (Fig. 2). 
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Fig. 4. Interaction of EK1 with HCoV. Side-by-side ribbon diagram and electrostatic surface representations illustrate that EK1 snugly fits into the hydrophobic groves 
formed between two adjacent HR1 helices of the 3HR1 core from MERS-CoV (A), SARS-CoV (B), and 229E (C). The EK1 peptide is shown as a green on white ribbon (left) 
and electrostatic (right) representations of corresponding HCoV 3HR1 cores. Hydrophobic surfaces are in whitish gray, basic in blue, and acidic in red. 


In all three structures, the EK1 peptide adopts a mixed secondary 
structure conformation (Fig. 4, A to C). The central region of EK1 
folds into a five-turn helix, which packs against two neighboring HR1 
helices via extensive hydrophobic interactions (3HRI cores of differ- 
ent HCoVs are illustrated as electrostatic surfaces, while residues on 
EK1 that are involved in hydrophobic packing are shown as stick 
models; residues L12, E15, M16, L19, A22, 123, L26, S29, and Y30 on 
EK1 locate within the five-turn helix of EK1 and form strong hydro- 
phobic interactions with the 3HR1 cores; Fig. 5A) and a few electro- 
static or polar interactions (residues E13, E15, K18, E20, E27, E28, 
and Y30 locate within the five-turn helix region of EK1 and interact 
with HRI residues through side chain-to-side chain hydrophilic 
interactions; Fig. 6A). The rest of EK1 adopts an extended conforma- 
tion, except that one extra turn is formed at the C-terminal end of 
EK1 (around E35) in the HR1(229E)-L6-EK1 structure (Fig. 6, right 
panel). In the extended region of EK1, polar interactions between side 
chain (HRI residues) and main chain (EK1 residues) dominate (com- 
pare Fig. 6, Band A). Most of these side chain-to-main chain polar 
interactions cluster at either end of the EK1 helical region, which in- 
terweave into extensive H-bond networks and likely help secure the 
EK1 helical region in the correct register (Fig. 6B). Hydrophobic res- 
idues in the extended region of EK1 also insert their bulky side 
chains into hydrophobic pockets on the surface of the 3HR1 cores 
(residues L2, I5, V7, L10, 131, and L36; Fig. 5A), further strengthen- 
ing the adhesion of the EK1 extended region onto the 3HRI cores. 

Both hydrophobic pockets and ridges exist on the surface of 3HR1 
cores (Fig. 5A). Correspondingly, we observed two kinds of hydro- 
phobic interactions between the EK1 and 3HRI cores. In particular, 
certain EK1 residues insert their hydrophobic side chains into 
pockets on the 3HR1 cores (Fig. 5A, shown as orange stick models 
and hereafter named as “burying residues”), and other EK1 residues 
pack their side chains against hydrophobic ridges on the 3HRI1 cores 
(Fig. 5A, shown as yellow stick models and hereafter named as 
“ridge-packing residues”). Note that the residues on HR1 that medi- 
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ate hydrophobic interactions with the burying and ridge-packing res- 
idues on EK] are conserved across all HCoVs [Fig. 5B (bottom panel), 
HRI residues shaded with orange and yellow], and most of the HR1 
residues involved in polar and electrostatic interactions with EK1 also 
manifest high levels of similarity across all HCoVs [Fig. 5B (bottom 
panel), HRI residues boxed in cyan and red]. Thus, EK1 would be 
able to form similar hydrophobic and hydrophilic interactions with 
HRIs from OC43, HKU1, and NL63, as we observed in the ho- 
mology models of HR1(OC43)-L6-EK1, HRI(HKU1)-L6-EK1, and 
HR1(NL63)-L6-EK1 (fig. $5). Together, extensive and highly con- 
served hydrophobic and hydrophilic interactions between EK1 and 
3HR1 cores endow EK1 with the ability to bind the 3HRI cores 
from different HCoVs and, hence, the capability of blocking the asso- 
ciation of different HR2s onto their corresponding 3HRI1 cores 
(Fig. 2). 

The EK1 peptide outcompetes self-derived HR2 of MERS-CoV in 
both cell-cell fusion and pseudovirus infection assays. We therefore 
also compared the interactions between EK1 and MERS-HR1 (Figs. 5A 
and 6, left panel) with those between MERS-HR2 and MERS-HR1 
(fig. S6A). Most hydrophobic interactions between EK1 and MERS- 
HRI are similar to those between MERS-HR2 and MERS-HR1 [com- 
pare Fig. 5A (left panel) and fig. S6A (left panel)]. Nevertheless, EK1 
formed more side chain-to-side chain hydrophilic interactions 
with MERS-HRI1 [compare Fig. 6A (left panel) and fig. S6A (middle 
panel)] than did MERS-HR2, which likely accounts for its better 
ICso in cell-cell fusion and pseudovirus infection assays (Fig. 2). 

Similar to HR1 from a-HCoVs, sequence alignment across 
HCoVs revealed that a 14-amino acid insertion also exists in the 
HR2 region of a-HCoVs (Fig. 1C). Such an insertion renders the 
HR2 helices of a-HCoVs four turns longer than those of B-HCoVs. 
The helical region of EK1 is only five turns long, much shorter 
than the nine-turn helix in 229E-HR2 (38). Nevertheless, the 
burying residues in the extended region of EK1 and their equivalents 
in 229E-HR2 all neatly insert their side chains into corresponding 
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Fig. 5. Hydrophobic packing between EK1 and 3HR1 hydrophobic cores are largely conserved across different EK1-HR1 complexes. (A) The 3HR1 cores are 
shown as electrostatic surfaces. At the positions where the hydrophobic surface on the 3HR1 core is deeply concave (pockets), EK1 residues that bury over 70% of 
their side-chain solvent accessible surface (SAS) into these pockets are shown as orange stick models. At the locations where the hydrophobic surface on the 3HR1 
core is relatively flat (ridges), EK1 residues that pack 50 to 70% of their side-chain SAS against these ridges are shown as light yellow stick models. EK1 residues are 
labeled, where those labeled in red form both hydrophobic and hydrophilic interactions with 3HR1 cores. (B) HR1 residues involved in interactions with EK1 are con- 
served across different HCoVs. EK1 and HR1 residues linked with dashed lines locate to the same layers on the 3HR1 triple helix. Burying EK1 residues are shaded 
orange, and ridge-packing EK1 residues are shaded light yellow. HR1 residues that mediate assembly of the 3HR1 cores are shaded orange, while those involved in 
ridge packing are shaded yellow. HR1 residues that mediate conserved side chain-to-side chain and side chain-to-main chain hydrophilic interactions with EK1 


residues are highlighted with cyan and red boxes, respectively. 


hydrophobic pockets (fig. S6C), although these residues are located 
in extended regions in EK1, while their equivalents are within 
o-helical regions in 229E-HR2 (fig. S6C). As a consequence, the 
hydrophobic interactions between EK1 and 229E-HRI1 are not 
compromised at all in comparison to those between 229E-HR2 
and 229E-HRI1 [compare Fig. 5A (right panel) and fig. S6B (left 
panel)]. Furthermore, a total of 12 polar and electrostatic in- 
teractions occur between EK1 and 229E-HRI1, including five 
side chain-to-side chain and seven side chain-to-main chain 
interactions [Fig. 6, A and B (right panel)]. Together, these ex- 
tensive hydrophobic and hydrophilic interactions preserve the 
high affinity of EK1 toward long HR1 HCoVs. That EK1 can form 
6-HB structures with both short (SARS-CoV and MERS-CoV) 
and long (229E) HR1s highlights its broad structural compatibility 
in accommodating HR1s from different HCoVs, thus consolidating 
its broad-spectrum inhibitory effect against pan-CoVs. 


DISCUSSION 

No effective and broad-spectrum anti-HCoVs drugs or vaccines are 
currently available in the clinic. Notwithstanding, the RBD of the 
S protein has been proposed as a promising target for the develop- 
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ment of specific antibodies and vaccines (39, 40). For example, the 
RBD-specific antibody CDC2-C2 exhibits inhibitory activity against 
MERS-CoV infection. Nevertheless, administration of the 
single-antibody CDC2-C2 could result in the emergence of es- 
cape mutations in MERS-CoV RBD (41). Unfortunately, the 
CoV S RBD is hypervariable throughout evolution, which has 
led to marked difference in host receptor usage in different 
HCoVs. Even when the same host receptor is used by different 
HCoVs, they frequently target different binding sites on the 
host receptor (42). Therefore, specific RBD-targeting antibodies 
or vaccines inevitably lack broad-spectrum activity against 
HCoV infection. For example, Menachery et al. (9) found that a 
SARS-CoV RBD-specific antibody could not protect mice from 
infection by the chimeric virus with the SHC014 S proteins, 
although SHC014 has high homology with SARS-CoV and can 
bind the same host receptor ACE2. The lag time between emerg- 
ing human CoV outbreaks and development of new prophylactic 
treatments or vaccines is of concern. Thus, there is an urgent need 
for the development of new, broad-spectrum drugs that 
target conserved sites in currently circulating and future emerging 
HCoVs so as to prepare for future outbreaks of yet unknown 
HCoVs. 
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Fig. 6. Highly conserved hydrophilic interactions between EK1 and 3HR1 cores. 
(A) Side chain-to-side chain hydrophilic interactions between EK1 and HR1. Side 
chains of EK1 and HR1 residues involved in these interactions are shown as cyan 
and gray stick models, respectively, and are similarly color-coded. At least four pairs 
of this hydrophilic interaction are conserved across different EK1-HR1 complexes: 
Y305' forms similar polar interactions with Q1009™""> (left), Q917°""> (middle), or 
T8177? (right); E27™' forms similar polar interactions with Q1009“"*® (left), Q917°""° 
(middle), or S8207"* (right); E15**! forms salt bridges with K1021™*"® (left), K929°""> 
(middle), or K83272** (right); and T8™“' makes similar interactions with N1028™="> 
(left), Q936°""> (middle), or Q83972”* (right). The conserved HR1 residues mentioned 
above are highlighted with cyan boxes. (B) Main chain-to-side chain hydrophilic 
interactions between EK1 and HR1.EK1 and HR1 residues involved in these interactions 
are shown as green and gray stick models, respectively. The main-chain atoms of 
HR1 residues are not shown for clarity. Similar to the side chain—-to-side chain inter- 
actions, main chain-to-side chain interactions are also highly conserved across 
different EK1-HR1 complexes. The HR1 residues involved in conserved interactions 
are highlighted with red boxes. 


In the early 1990s, a series of peptides derived from the HIV-1 
gp41 HR2 (or CHR) domain, such as SJ-2176 (19), DP-178 (later 
T20) (20), and C34 (21), were reported to have highly potent inhib- 
itory activity against HIV-1 gp41-mediated membrane fusion and 
HIV-1 infection with ICs 9 values at low nanomolar levels. Subse- 
quently, numerous virus fusion inhibitory peptides overlapping the 
HR2 sequences of class I membrane fusion proteins from other 
enveloped viruses, including RSV (22), Ebola virus (23), paramyxo- 
viruses SV5 (24), and Nipah virus (25), were also reported. In 2003, 
Bosch et al. (26) discovered that a soluble HR2 peptide derived from 
the MHV S protein HR2 region inhibited virus cell entry, suggesting 
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that the S glycoprotein of CoV is a class I fusogen with the ability to 
form 6-HB during the S protein-mediated membrane fusion pro- 
cess. In 2004, the Bosch (43) and Jiang (18) groups independently 
reported that peptides derived from the HR2 region of SARS-CoV S 
protein could interact with the peptides from its HR1 region to form 
6-HB and inhibited SARS-CoV S protein-mediated membrane 
fusion and SARS-CoV infection with moderate potency. In 2014, 
Lu et al. (27) reported that peptide derived from the MERS-CoV 
S-HR2 could competitively inhibit 6-HB formation, thereby prevent- 
ing fusion of the virus with host cells. 

In our previous studies, we found that the MERS-CoV-specific 
fusion inhibitor MERS-HR2P could not block SARS-CoV pseudo- 
virus infection and did not display any broad-spectrum inhibitory 
activity. Thus, it was unclear whether 6-HB was a good target site 
for the development of a broad-spectrum anti-HCoV inhibitor. To 
identify a pan-HCoV fusion inhibition target site and preferably also 
to develop a potential pan-HCoV inhibitor against infection of 
multiple HCoVs in the human respiratory tract, we successfully 
established multiple HCoV S-mediated cell-cell fusion assays to de- 
termine the cross-inhibitory spectrum between HR1Ps and HR2Ps, 
which are derived from the HR1 and HR2, respectively, of HCoVs. 
Unexpectedly, we found that OC43-HR2P harbored broad-spectrum 
binding activity to all HR1Ps and manifested fusion inhibitory 
activity against all HCoV S-mediated cell-cell fusion. We then further 
optimized the OC43-HR2P to improve its antiviral activity and 
solubility. The optimized peptide EK1 exhibited broad and potent 
fusion inhibitory activity against infection by HCoVs and SL-CoVs 
on various in vitro models, i.e., cell-cell fusion assays and pseudo- 
typed or live viral infection assays. We assessed the prophylactic and 
protective effects of EK1 through intranasal administration on OC43 
and MERS-CoV mouse models. EK1 exhibited effective preventive 
and protective effects in these mouse models with an acceptable in 
vivo safety profile. Furthermore, the in vivo pharmacokinetic pro- 
filing and safety studies helped us acquire useful data for potential 
human clinical trials in the future. 

The potency of HR2-derived peptides in inhibiting corresponding 
virus-cell fusion varies significantly. For example, anti-HIV-1 peptide 
T20 is about 900-fold more potent than the anti-SARS-CoV pep- 
tide SC-1, yet only 30-fold more potent than the anti-MERS-CoV 
peptide HR2P (18, 27, 44). It is known that HIV-1 and SARS-CoV 
each enter their target cells mainly through plasma and endosomal 
membrane fusion, respectively, while MERS-CoV could infect its 
target cells via both plasma and endosomal membrane fusion 
(18, 27, 44). The vast differences seen in potency thus suggest that 
only a limited number of the EK1 peptides can get into the endo- 
some to inhibit virus-cell fusion. Therefore, increasing the cell per- 
meability of the fusion inhibitory peptides is expected to enhance 
the potency of these peptides. We have recently reported that addi- 
tion of hydrocarbon stapling or a palmitic acid group (C16) to the 
MERS-CoV fusion inhibitory peptides significantly improves 
their antiviral potency and pharmacokinetic properties (45, 46). 
In the near future, we will use the similar approaches to improving 
the antiviral activity of EK1. 

Recent studies have reported that several emerging SL-CoVs, 
including WIV1, SHCO014, and Rs3367, have the potential to infect 
humans with high virulence by further evolution or direct gene 
recombination with SARS-CoV (9, 11). We found that the HR se- 
quences of these SL-CoVs are the same as those of SARS-CoV, and 
we predicted that EK1 could prevent infection by these SL-CoVs. As 
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expected, EK1 efficiently inhibited SL-CoV S protein-mediated cell- 
cell fusion. Consistent with previous studies (9, 11), we herein also 
failed to establish a SHC014 pseudovirus using an HIV-1 backbone 
vector; therefore, we assessed the antiviral effect of EK1 only on 
WIV1 and Rs3367 pseudovirus infection and found that EK1 ex- 
hibited effective inhibitory activity against these two SL-CoVs. In 
addition, some animal CoVs within the CoV family are very close to 
HCoV, such as bovine coronavirus (BCoV) and MHV. The HR se- 
quences of these CoVs are very similar to those of OC43, indicating 
that peptide EK1 will likely also be effective on these viruses. 

One intriguing question is why only EK1 and the HR2 from 
OC43 manifested “pan-CoV” inhibitory activity, while HR2 from 
other HCoVs did not. HR2s from a-HCoVs are substantially longer 
than those from B-HCoVs (Fig. 1C), and docking of HR2s from 
229E and NL63 onto the 3HR1 cores of B- HCoVs resulted in severe 
steric clashes (figure not shown). It is thus understandable why the 
HR2s from 229E and NL63 would not have broad-spectrum activity 
toward B-HCoVs (Fig. 1D). 

At the same time, sequence alignment between EK1 and HR2s 
from OC43, SARS-CoV, and MERS-CoV provided clues on the very 
limited inhibition breadth of SARS-HR2 and MERS-HR2. Although 
most residues that mediated hydrophobic (shaded orange and yellow) 
and side chain-to-side chain hydrophilic interactions (boxed in cyan) 
are conserved between EK1 and different HR2s, sequence variability 
does exist at positions V7, Y30, and some other ridge-packing resi- 
dues (fig. S7A). We thus focused our analysis on these residues and 
their surrounding environments in different HCoVs 3HRI cores. In 
both solved structures of MERS(HR1)-EK1 and SARS(HR1)-EK1, 
V7'"' fits its side chain neatly into the hydrophobic pockets on these 
two 3HRI cores (fig. S7B, second and fourth panels). Nevertheless, 
in the previously reported structure of the SARS-CoV fusion core, 
A1156°48S, the equivalent of V7§*! in SARS-CoV, would occupy the 
same cavity loosely (fig. S7B, third panel), while T1257M=®S, the 
equivalent of V7=*! in MERS-CoV, could not fully bury its side 
chain owing to the polar nature of the Thr residue (fig. S7B, first 
panel). Thus, Val, in EK1 and OC43-HR2, is more effective in fitting 
the hydrophobic cavities at this position than Ala (seen in SARS- 
HR2) or Thr (seen in MERS-HR2) (fig. S7B). 

In total, we observed four ridge-packing interactions between the 
EK1 and 3HRI cores. The ridge to which $29Y30""' binds is sur- 
rounded by polar and positively charged residues in MERS-CoV and 
SARS-CoV but surrounded by polar and negatively charged residues 
in 229E [fig. S7C (top panel), side chains of surrounding residues 
are shown as stick models]. At this position, the polar pair “Ser-Tyr” 
in EK1 and OC43-HR2 appears to be more suitable than the “Ser-Leu” 
pair from SARS-HR2 in accommodating the opposite electrostat- 
ic environments in different HCoVs. Moreover, the aromatic ring 
of Tyr could also form additional polar interactions with HR1 side 
chains in different HCoVs (Fig. 6A). Likewise, the ridges to which 
E15M16"" binds also exhibited better shape complementarity to 
Met (in EK1 and OC43-HR2) than to Ile (in SARS-HR2). Collec- 
tively, it seems that EK1 and OC43-HR2 have more optimal amino 
acid choices at all the abovementioned critical positions, which likely 
accounts for their unique pan-CoV inhibitory activity. 

In its natural state, the S protein that is present on the CoV surface 
is inactive. After receptor binding target cell proteases activate the S 
protein by cleaving the exposed enzyme target sites, leaving the S2 
subunit free to mediate viral fusion and entry. CoV can enter the 
target cell via two pathways: one is the endocytosis pathway and the 
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other is direct fusion on the cellular surface. For example, when 
SARS-CoV enters the target cell through the endocytosis pathway, its 
S protein can be cleaved and activated by the pH-dependent cysteine 
protease cathepsin L in the endosome. On the other hand, recent 
studies have consistently reported that the SARS-CoV S protein can 
also be cleaved and activated by transmembrane protease serine 2 
(TMPRSS2) and human airway trypsin-like protease, which are 
located on the cell surface, thus activating and allowing the S protein 
to mediate SARS-CoV infection at the plasma membrane (9-11). In 
accord with such a finding, Matsuyama et al. (47) reported that the 
HR2 peptide efficiently inhibited SARS-CoV entry into cells, while 
lysosome-tropic reagents failed to inhibit at all. Similarly, TMPRSS2 
also has the capacity to promote the entry of MERS-CoV through 
bypassing the endocytosis pathway and directly activating its S pro- 
tein on the cellular surface (48). Recently, several studies have reported 
that TMPRSS2 is highly expressed on human respiratory epithelial 
cells surface and was even associated with several CoV receptors, 
such as ACE2 and DPP4 (48-50). Hence, the plasma membrane 
fusion pathway seems a preferred choice for HCoV infection in the 
human respiratory tract. Consistently, our previous study found that 
MERS-CoV-specific fusion inhibitor HR2PM2 effectively inhibited 
MERS-CoV infection in vivo by intranasal administration (33). Other 
studies have also reported that the current clinical isolates of 229E 
and OC43 are very sensitive to cell surface TMRRSS2 but not to endo- 
somal cathepsins (51, 52). Overall, for current circulating HCoVs or 
emerging HCoVs, the cell surface fusion pathway in human respi- 
ratory tract would appear to be very important. Therefore, peptide 
fusion inhibitors and the strategy of intranasal administration are 
excellent choices for preventing HCoV infection via the airway, 
which is a key site for HCoV to rapidly establish infection and widely 
spread to other organs. 

Currently, circulating HCoVs pose a potential threat to humans; 
moreover, it is almost certain that other zoonotic CoVs will be trans- 
mitted to humans in the future. The availability of HCoV-specific 
drugs with broad-spectrum inhibitory activity is therefore important 
for the prevention and control of a future HCoV epidemic. The EK1 
fusion inhibitor peptide targeting the conserved site in the spike HR1, 
but not the hypervariable RBD region, has potent and broad inhibi- 
tory activity against multiple HCoV infections. EK1 through nasal 
administration exhibited effective and broadly anti- HCoV activity 
with a satisfactory safety profile in vivo. Hence, EK1 is a promising 
candidate for further development as an antiviral agent against 
infection of multiple HCoVs, especially for use in infants and the 
elderly, as well as immunocompromised patients, who would be more 
vulnerable to HCoV infections (53-56). Meanwhile, this study pro- 
vides clues and methods for the development of peptide fusion inhibi- 
tors with potency and breadth in inhibiting infections by other 
highly pathogenic enveloped viruses with class I membrane fusion 
proteins, such as Ebola and Marburg viruses, Hendra and Nipah 
viruses, and influenza viruses. 


MATERIALS AND METHODS 

Cells, viruses, and peptides 

Cell lines Huh-7, Vero E6, HCT-8, Calu-3, A549, LLC-MK2, 293 T, 
and 293 T/ACE2 were obtained from American Type Culture Col- 
lection (ATCC; Manassas, VA, USA). All cell lines were cultured in 
Dulbecco’s modified Eagle’s medium (DMEM) with 10% fetal 
bovine serum (FBS). The HCoV-OC43 (VR-1558) and HCoV-229E 
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(VR-740) strains were obtained from the ATCC. MERS-CoV- 
EMC/2012 was originally provided by H. Feldmann (National 
Institutes of Health, National Institute of Allergy and Infectious 
Diseases, Rocky Mountain Laboratories, Hamilton, MT) and R. A. 
Fouchier (Erasmus Medical Center, Rotterdam, Netherlands). 
PMD2.G-VSV-G, pcDNA3.1-MERS-S, pcDNA3.1-SARS-S, pcDNA3.1- 
WIVL-S, pcDNA3.1-Rs3367-S, pCDNA3.1-SHC014-S, pNL4-3.Luc.R-E, 
and pAAV-IRES-EGFP plasmids and the NL63 strain (Amsterdam) 
were preserved in our laboratory. The DNA sequence of the S pro- 
tein of OC43 with a deletion of 17 amino acids in the C terminus 
was synthesized. DNA encoding 229E-S or NL63-S-18 (containing a 
C-terminal 18-amino acid deletion) was provided by F. Li. All peptides 
were synthesized by KareBay Biochem with a purity of >95% (tested 
by high-performance liquid chromatography). The EK1-scrambled se- 
quence was “LK VLLYEEFKLLESLIMEILEYQKDSDIKENAEDTK.” 


CoV phylogenetic analysis 

Phylogenetic trees were constructed using the MEGA6.06. Acces- 
sion numbers used for phylogenetic analysis are as follows: BCoV 
(KX982264.1), PEDV (NC_003436), HCoV-229E (AAK32191.1), 
HCoV-NL63 (AVA26873.1), MHV (AC_000192.1), IBV (KY421672), 
HCoV-OC43 (CAA83661.1), HCoV-HKU1 (DQ415904), BtCoV-HKU4 
(NC_009019.1), BtCoV-HKU5 (NC_009020), SW1 (EU111742.1), 
MERS-CoV (AID55097.1), SARS-CoV (ABD72979.1), BtSCoV-SHC014 
(KC881005), BtSCoV-Rs3367 (KC881006), BtSCoV-WIV 1 (KF367457), 
BtCoV-HKU9 (EF065516), HKU11 (FJ376620), HKU12 (NC_011549.1), 
HKU13 (NC_011550.1), and PDCoV (KX022605.1). 


Expression and purification of fusion protein HR1-L6-EK1 
The coding sequence of EK1 peptide was individually fused to the 3’ 
end of the HR1 domain from SARS-CoV, MERS-CoV, and HCoV- 
229E (residues 892 to 970, 984 to 1062, and 785 to 873, respectively) 
through a six-amino acid linker (L6: SGGRGG). The resulting 
sequences encoding different HR1-L6-EK1 fusion proteins were 
then subcloned into a modified pET-28a vector, which contains a 
Hisg-SUMO tag upstream of the multiple cloning site. The resulting 
constructs, pET-28a- Hisge-SUMO-HR1-L6-EK1, were then expressed 
in Escherichia coli BL21 (DE3) at 16°C overnight in LB medium. We 
initially purified these fusion proteins using His-Talon resin (Clon- 
tech). Eluted fractions from the His-Talon column were then mixed 
with Ulpl enzyme [1:100 (w/w)] and dialyzed against buffer A 
[20 mM tris-HCl (pH 8.0), 150 mM NaCl, and 1 mM dithiothreitol] 
at 4°C overnight. After SUMO tag cleavage, the samples were re- 
loaded onto the His-Talon column, and flow-through fractions 
containing untagged HR1-L6-EK1 were pooled, concentrated, 
and gel-filtered in buffer B [20 mM tris (pH 8.0) and 150 mM 
NaCl] on a HiLoad 16/60 Superdex 75 (GE Healthcare) column. 
Peak fractions that contain HR1-L6-EK1 trimer were pooled and con- 
centrated to 10 mg/ml through centrifugation (EMD Millipore). 


Crystallization and structure determination 

Initial crystallization conditions were screened using the sitting- 
drop vapor diffusion method on a Mosquito crystallization robot 
(TTP Labtech) at 20°C. Diffraction-quality crystals of the HR1(SARS)- 
L6-EK1 fusion protein were obtained using the hanging-drop vapor 
diffusion method by mixing equal volumes (0.5 ul) of protein solu- 
tion (5 mg/ml) and reservoir solution [0.2 M Li,SO,, 0.1 M bis-tris 
(pH 5.5), and 25% PEG3350]. Diffraction-quality crystals of the 
HR1(MERS)-L6-EK1 and HR1(229E)-L6-EK1 fusion proteins were 
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obtained through similar methods from 0.2 M KSCN (pH 7.0), 20% 
PEG3350, and 0.05 M MgCl, and 0.1 M HEPES (pH 7.5) and 30% 
PEG550MME, respectively. Diffraction data were collected at beam- 
line BL19U1 of Shanghai Synchrotron Radiation Facility (SSRF), 
China, and processed with the HKL3000 program (57). A summary 
of the data collection statistics is provided in table $3. The structure 
of HR1(SARS)-L6-EK1 was solved by molecular replacement, as im- 
plemented in the PHASER program of PHENIX (58). The programs 
used a SARS-CoV fusion core structure [Protein Data Bank (PDB): 
2bez] as the search model. For structure determination of HR1 
(MERS)-L6-EK1 and HR1(229E)-L6-EK1, the MERS-CoV fusion core 
[PDB: 4mod] and SARS-CoV fusion core [PDB: 2bez] were used as 
search models. The structural models were further improved by cy- 
cles of manual building and refinement using the COOT (59) and 
PHENIX (58) programs. The quality of these models were analyzed 
with MolProbity (60). A summary of the structure refinement sta- 
tistics is also given in table $3. The figures were all prepared using 
the PYMOL program (The PyMOL Molecular Graphics System, 
Version 2.1, Schrédinger LLC). The electrostatic calculations were 
performed with PDB2PQR (61). 


Homology model building 

Models for HR1(OC43)-L6-EK1, HR1(HKU1)-L6-EK1, and 
HR1(NL63)-L6-EK1 were derived by homology modeling using 
Swiss Model website (62). The template for HR1(OC43)-L6-EK1 
and HR1(HKU1)-L6-EK1 were obtained from the crystal structures 
of HR1(SARS)-L6-EK1 and HR1(MERS)-L6-EK1, and the template 
for HR1(NL63)-L6-EK1 proteins was obtained from the crystal 
structure of HR1(229E)-L6-EK1. To relax and stabilize the interac- 
tion between EK1 and corresponding HRI, initial models were 
optimized by performing energy minimization, followed by a 5-ns 
molecular dynamics simulation using Schrédinger Suite 2017-4 
(www.schrodinger.com). The simulation systems were solvated with 
full-atom TIP3P water, containing Cl’ and Na’ ions at a concentra- 
tion of 0.15 M to mimic physiological ionic strength. During the 
simulation, temperature T and pressure P were kept constant, at 
310 K and 1 atm, respectively. 


Circular dichroism spectroscopic analysis 

Circular dichroism spectra (195 to 260 nm) were collected on a 
J-815 spectropolarimeter (JASCO Inc., Japan) to evaluate the second- 
ary structure of the peptides and their complexes (63, 64), individual 
peptides or complexes dissolved in PBS with the final concentration 
at 10 uM. Thermal denaturation of peptide complexes was monitored 
from 20° to 100°C at 222 nm with a thermal gradient of 5°C/min. 
The midpoint of the T,, values was acquired by JASCO software util- 
ities. 


Inhibition of pseudotyped virus infection 

A pseudovirus bearing CoV S protein or VSV-G protein and a 
defective HIV-1 genome that expresses luciferase as reporter was 
produced in 293 T cells, as previously described (27), and its titer was 
quantitated by using HIV-1 p24 ELISA. The pseudovirus was then 
used to infect target Huh-7 cells (or ACE2/293 T cells for pseudo- 
typed SARS-CoV) (10 * per well in 96-well plates) in the presence or 
absence of the test peptide at the indicated concentration. Twelve 
hours after infection, culture medium was refreshed and then incu- 
bated for an additional 48 hours, followed by washing cells with 
PBS, lysing cells with lysis reagent (Promega), and transferring 


11 of 15 


6LOZ ‘OL |udy uo /Bio bewseousi0s ssoueApe//:dy1y Wos) papeojuMOG 


SCIENCE ADVANCES | RESEARCH ARTICLE 


the cell lysates to 96-well Costar flat-bottom luminometer plates 
(Corning Costar) for the detection of relative light units using 
the Firefly Luciferase Assay Kit (Promega) and an Ultra 384 lumi- 
nometer (Tecan). 


Establishment of multiple HCoV S-mediated cell-cell fusion 
As previously described (27), 293 T effector cells were transfected 
with plasmid pAAV-IRES-EGFP encoding the EGFP (293 T/EGFP 
cells) or plasmid pAAV-IRES-S-EGFP encoding the corresponding 
HCoV S protein (293 T/HCoV S/GFP cells) as the effector cells. 
Huh-7 cells, expressing various HCoV receptors on the mem- 
brane surface, were used as target cells, as described below. 

1) MERS-CoV S-mediated cell-cell fusion: Effector cells (293 T/ 
MERS-CoV/GEFP) and target cells (Huh-7 cells) were cocultured in 
DMEM containing 10% FBS, at 37°C for 2 hours; 

2) 229E S-mediated cell-cell fusion: Effector cells (293 T/229E/ 
GFP) and target cells (Huh-7 cells) were cocultured in DMEM con- 
taining 10% FBS, at 37°C for 4 hours; 

3) SARS-CoV and SL-CoV S-mediated cell-cell fusion: Effector 
cells (293 T/SARS-CoV/GFP or 293 T/SL-CoV/GFP) and target cells 
(Huh-7 cells) were cocultured in the presence of trypsin (80 ng/ml) 
in DMEM without FBS, at 37°C for 4 hours; 

4) OC43 or NL63 S-mediated cell-cell fusion: Effector cells (293 T/ 
HCoV-OC43/GFP or 293 T/HCoV-NL63/GFP) and target cells 
(Huh-7 cells) were cocultured in the presence of trypsin (80 ng/ml) 
in DMEM without FBS, at 37°C for 4 hours. 

Five fields in each well were randomly selected for counting the 
fused and unfused cells. The fused cells are at least twice as large as 
the unfused cells, and the fluorescence intensity in the fused cell 
became weak as a result of the diffusion of enhanced green fluores- 
cent protein (EGFP) from one effector cell to target cells (see figs. 
S1D and S2C). The percentage of cell-cell fusion [(number of the 
fused cells/number of the fused and unfused cells) x 100%] was then 
calculated. 


Inhibition of HCoV S-mediated cell-cell fusion 

The inhibitory activity of a test peptide on HCoV S-mediated cell- 
cell fusion was determined, as previously described (27). Briefly, 
effector cells (293 T/S/GFP) and target cells (Huh-7 cells) were co- 
cultured in the presence or absence of a test peptide at the indicated 
concentrations for fusion. After counting the fused and unfused 
cells, the percentage of cell-cell fusion was calculated, as described 
above. The percent inhibition of cell-cell fusion was calculated 
using the following formula as described elsewhere (27): [1 — (E- N)/ 
(P — N)] x 100%. “E” represents the percentage of cell-cell fusion in 
the experimental group. “P” represents the percentage of cell-cell 
fusion in the positive control group, where 293 T/HCoV S/EGFP 
cells were used as effector cells, to which only PBS was added. “N” 
is the percentage of cell-cell fusion in negative control group, 
where 293 T/EGFP cells were used as effector cells. The ICs9 was 
calculated using the CalcuSyn software provided by T. C. Chou 
(65). Samples were tested in triplicate, and all those experiments were 
repeated twice. 


HCoV virion-based fusion assay 

An HCoV virion-based fusion assay was performed, as described 
elsewhere (34, 35). Briefly, Huh-7 cells and ACE2-293 T cells were 
used as target cells for the entry of MERS-CoV Blam-Vpr virions 
and SARS-CoV Blam-Vpr virions, respectively. These target cells 
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were cultured at 37°C for 5 hours in six-well plates (with virions 
containing BlaM-Vpr equivalent to 80 ng of p24-Gag per well) in 
the presence or absence of EK1 at the indicated concentrations. The 
cells were washed with PBS, resuspended in 500 ul of DMEM, and 
then incubated with CCF4-AM substrate at room temperature for 
2 hours, as described by the manufacturer (Invitrogen, Germany). 
Last, the cells were monitored via flow cytometry. After the virion 
fusion with the target cell, the CCF4-AM (emission at 520 nm) 
substrate could be cleaved by BlaM-Vpr into CCF4 (emission at 
447 nm). Flow cytometric data were collected with a BD FACSDIVA 
and analyzed with FlowJo. 


Inhibition of live HCoV’s replication 

The inhibitory activity of peptides against OC43 replication in 
HCT-8 cells was assessed, as described elsewhere (66). Briefly, 
100 TCIDs9 of OC43 was mixed with a test peptide at graded con- 
centration and incubated at 37°C for 30 min. The mixture was then 
applied in triplicate onto the monolayer of HCT-8 cells grown in a 
96-well microtiter plate. On day 5 after infection, viral titer in the 
culture medium was tested, and TCIDs9 was calculated on the basis 
of the cytopathic effect (CPE) (67). The inhibitory activity of the 
tested peptide against 229E replication in A549 cells and NL63 rep- 
lication in LLC-MKz2 cells was evaluated in a similar way, as de- 
scribed above. 

The inhibitory activity of peptide against MERS-CoV replication 
was tested in Calu-3 cells using a modified standard microneutraliza- 
tion assay, as previously described (68). Briefly, 60 ul of a serially 
twofold diluted peptide was incubated with 60 wl (120 TCIDs 9) of 
MERS-CoV in MEM medium supplemented with 2% FBS (M-2 
medium) in duplicate wells of 96-well plates for ~60 min at room 
temperature. One hundred microliters of the peptide/MERS- 
CoV mixtures was then transferred into confluent Calu-3 cells 
grown in 96-well plates. Wells of Calu-3 cells cultured with M-2 
medium with and without virus were included in these assays as pos- 
itive and negative controls, respectively. While the advanced CPE of 
Calu-3 cells could develop within 24 to 36 hours in response to 
MERS-CoV infection at higher multiplicities of infection (MOIs) of 
1 or 0.1 (69), we did not observe any prominent formation of CPE 
until ~ 60 to 72 hours after infection at an estimated MOI of ~ 0.001. 
Hence, to more accurately measure the efficacy of peptide inhibitors 
against MERS-CoV infection, we harvested the supernatants at 72 hours 
and quantified infectious virus titers by the standard Vero E6-based 
infectivity assays and expressed the titers as logj9 TCIDs59/ml. 


Mouse infection studies 
Pregnant Balb/c mice (18 days) purchased from the Department 
of Laboratory Animal Science, Fudan University were separated 
into four groups after delivery of their offspring. Eleven newborn mice 
were chosen for each group. Mice in the prevention and treatment groups 
were intranasally administered peptide (5 mg/kg in 2 ul of PBS) 30 min 
before or after intranasal challenge with a viral dose of 107 TCIDso 
(in 2 ul DMEM). Mice in the viral control group and the normal control 
group were intranasally administered with 2 ul of PBS 30 min before 
viral challenge or without viral challenge. Mouse survival rate and 
body weight variations were recorded up to 2 weeks after infection (70). 
On day 5 after infection, five mice in each group were randomly 
selected for euthanasia to collect and assess the viral titer in mouse brain. 
Briefly, mice expressing human DPP4 were randomly separated 
into three groups (n = 6): viral control, prophylactic group, and therapeutic 
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group. In the viral control and therapeutic groups, each mouse was 
challenged intranasally with 10* TCIDsp [i-e., 1000 LDsp (median le- 
thal dose)] of MERS-CoV in a volume of 60 ul, 30 min after viral 
challenge, and then treated via the intranasal route with 200 ug of 
EK1 (in 50 wl of PBS) or 50 ul of PBS, respectively. Prophylactic group 
mice were pretreated via the intranasal route with 200 ug of EK1 per 
mouse 30 min before viral challenge. Challenged mice were moni- 
tored daily for weight loss and mortality. In the interim, two mice 
were randomly euthanized (two of six) in each group at day 2 after 
challenge to determine the lung virus titers by CPE on Vero E6 and 
expressed as TCIDs0/ml. 


In vivo fluorescence imaging 

Six Balb/c female mice (8 weeks old) were randomly assigned to two 
groups and nasally administered with 20 ug of Cy5-EK1 in 50 wl of 
PBS (n = 3) or 50 ul of PBS without peptide (n = 3) as control for 
background fluorescence measurement. One hour later, mice were 
imaged for distribution of Cy5-EK1 by the IVIS Lumina K Series I] 
in vivo imaging system (PerkinElmer, Waltham, MA, USA). Mice 
were euthanized by intraperitoneal injection of sodium pentobarbital 
to obtain the lungs, livers, kidneys, spleens, and hearts for imaging 
and measuring their fluorescence value. The relevant radiant effi- 
ciency (P s! cm” sr!) (uw cm’) was then calculated by Living 
Image 4.4 software (56, 71, 72). 


Cytotoxicity assay 

Cytotoxicity of the peptides to the cells (293 T, 293 T/ACE2, Huh-7, 
A549, LLC-MK2, and Calu-3 cells) was tested by using the Cell 
Counting Kit-8 (CCK-8; Dojindo, Kumamoto, Japan). Briefly, each 
cell type was seeded into the wells of a 96-well microtiter plate (10 * 
per well) and incubated at 37°C for 12 hours, replacing medium 
with DMED containing EK1 at graded concentrations. After incu- 
bation at 37°C for 2 days, CCK-8 solution (10 ul per well) was 
added, followed by an additional incubation for 4 hours. The ab- 
sorbance was measured at 450 nm. 


In vivo safety of EK1 through intranasal administration 

As previously described (27), female Balb/c mice (8 weeks old) were 
assigned randomly to three groups (m = 5) and continuously ad- 
ministered with PBS, low-dose EK1 (20 mg/kg), or high-dose EK1 
(100 mg/kg) by intranasal route every day for 1 week. The body 
weight changes were monitored for the following 2 weeks. ALT and 
creatinine in the sera of each group of mice were measured using 
the ALT and creatinine assay kits (NJJCBIO, Nanjing, China) before 
the first treatment and 1, 3, and 5 days after the final treatment. The 
titer of immunoglobulin G sera to EK1 in each group of mice was 
evaluated by ELISA 2 weeks after the final treatment. Four weeks 
after administration, mice in each group mice were euthanized to 
harvest the lungs, livers, kidneys, and spleens for hematoxylin and 
eosin staining. The animal studies were approved by the Institutional 
Laboratory Animal Care and Use Committee at Fudan University 
(20180302-019). 


Statistical analyses 

Analyses of independent data were performed by Student’s unpaired 
two-tailed ¢ test. Statistical analyses were carried out using GraphPad 
Prism 6.0. P values less than 0.05 were considered significant; *P < 0.05, 
**P < 0.01, ***P < 0.001, and ****P < 0.0001. The ICs9 was calculated 
using the CalcuSyn software (65). 
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SUPPLEMENTARY MATERIALS 
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